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We study the consequences of in-plane (D p ) and out-of-plane (D z ) Dzyaloshinsky-Moriya (DM) 
interactions on the thermodynamic properties of spin-i Heisenberg model on the kagome lattice 
using numerical linked cluster expansions and exact diagonalization, and contrast them with those 
of other perturbations such as exchange anisotropy and dilution. We find that different combinations 
of the DM anisotropies lead to a wide variety of thermodynamic behavior, which are quite distinct 
from those of most other perturbations. We argue that the sudden upturn seen experimentally in 
the susceptibility of the material ZnCu3(OH)6Cl2 can be understood in terms of Dzyaloshinsky- 
Moriya anisotropies with D P >\D Z \. We also show that the measured specific heat of the material 
puts further constraints on the allowed DM parameters. 

PACS numbers: 75.10. Jm,05.50.+q,05.70.-a 
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I. INTRODUCTION 

The possible realization of exotic states of matter has 
always been at the forefront of research interest in con- 
densed matter physics. One class of states, which has 
received considerable interest over the last few decades, 
consists of quantum spin liquids^ In these states, no mag- 
netic order or other symmetry breaking occurs as the 
temperature is lowered, while the system may or may 
not exhibit a gap in its excitation spectra. Two variants 
of the quantum spin liquid have received particular at- 
tention recently: (i) a topological spin liquid, which has 
a spin gap and a topological order^ and (ii) an alge- 
braic spin liquid, where there is no spin gap and spin- 
spin correlations decay as a power law.— & Kagome lattice 
antiferromagnets are potential candidates for both these 
types of spin liquids. 

Recently, the newly synthesized herbertsmithite 6 
ZnCu3(OH)eCl2 has brought tremendous excitement to 
the field. For this rare mineral, in which the spin-^ cop- 
per atoms form a kagome lattice, like the one depicted 
in Fig. [I] no magnetic order is observed down to 50 mK 
(Refs.0H!|l3) even though the exchange constant is ap- 
proximately 170 K^i 

The spin-i kagome lattice Heisenberg model (KLHM) 
has been extensively studied using series expan- 
sions and exact diagonalization of finite-size periodic 
clusters-Hil&IMM&iLI&iS Exact diagonalization (ED) 
studies suggest that this model does not exhibit any mag- 
netic order but possibly has a small spin gap ~ J/20. In 
addition, in finite systems, it has been found that there 
are a large number of singlet states below the spin gap 
and that their number grows with the system size,— ^ 
which indicates that in the thermodynamic limit, non- 
magnetic excitations may develop a continuum beside the 
ground state. 

Considering the above theoretical results, the experi- 
mental behavior of the recently synthesized kagome sys- 



tems, ZnCu3(OH) 6 Cl2, is highly unexpected^ At high 
temperatures, the inverse susceptibility data was found 
to obey a Curie- Weiss law, with an effective Curie- Weiss 
constant of about 300 K. However, no spin gap was seen 
either in the susceptibility, the specific heat, the neu- 
tron spectra, or in the nuclear spin-lattice relaxation T\ , 
down to temperatures below 100 mK. In addition, at the 
lowest temperatures, the susceptibility saturates to very 
high values and the specific heat shows power-law behav- 
ior in temperature. The latter is suppressed by magnetic 
fields, showing it to be magnetic in origin. 




FIG. 1: (Color online) The kagome lattice. The embedded 
finite clusters are used (with periodic boundary conditions) 
in the ED study of DM interactions. 

A possible interpretation of the substantial rise in the 
susceptibility seen experimentally at low temperatures is 
that it is due to impurity spins outside the kagome planes, 
possibly caused by substitutions of non magnetic Zn sites 
with Cu. 21 This idea has been reinforced in recent numer- 
ical analysis of Misguich and Sindzingre^ They find that 
the deviation of the experimental data from the kagome 
lattice Heisenberg model at least down to T«40 K can be 
understood in terms a small concentration of impurities 
provided the impurity contribution also has a weak ferro- 
magnetic Curie- Weiss constant. We, on the other hand, 
have argued that the sharp rise in the susceptibility seen 
experimentally is intrinsic and it is related to the presence 
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of Dzyaloshinsky-Moriya (DM) interactions. 11 The latter 
is consistent with the experimental observation that the 
muon shift K tracks the bulk susceptibility %. 

The two different explanations can be distinguished 
by studying single crystals, whereby the susceptibilities 
along different crystallographic axes can be investigated. 
The explanations based on DM interactions lead to sub- 
stantial temperature dependent anisotropy in the sus- 
ceptibility, whereas the impurity contributions should be 
isotropic. In the absence of single crystals, recent NMR 
work of Imai et al^ may already provide some resolu- 
tion of the issue. Imai et al. find that the NMR spec- 
tra progressively broaden as one goes to lower tempera- 
tures. While the temperature dependence of the median 
of the broadened spectra resembles the sharp upturn seen 
in bulk susceptibility measurements, the spectra at the 
edges do not show this upturn. Rather, they show a 
saturation and eventually a downturn with lowering of 
temperature as expected in antiferromagnets, when short 
range order sets in. Indeed, we will see here that for cer- 
tain choices of the DM parameters, susceptibilities along 
different crystallographic axes show the two different be- 
haviors observed by Imai et al. Thus, one interpreta- 
tion of the NMR experiments is that they are observing 
the susceptibility along different axes, due to the differ- 
ent alignment of the powders with respect to the applied 
field. 

In this work, we further develop the DM theory by 
presenting a detailed study of the effects of in-plane 
(Dp) and out-of-plane (D z ) Dzyaloshinsky-Moriya inter- 
actions, as well as other perturbations such as easy-plane 
and easy-axis exchange anisotropies and quenched dilu- 
tion, on the thermodynamic properties of the KLHM. 
We use ED and the triangle-based numerical linked clus- 
ter method (NLC)22 to compute the uniform suscepti- 
bility, entropy, and specific heat. We find that unlike 
exchange anisotropy or dilution, the effects of DM inter- 
actions on the susceptibility can set in quite abruptly 
as a function of temperature. Furthermore, different 
choices of the DM parameters can lead to a wide range 
of behavior for the susceptibility. In particular, the 
abrupt upturn in the susceptibility seen experimentally in 
ZnCu3(OH)6Cl2 around 75 K can be understood in terms 
of Dzyaloshinsky-Moriya anisotropies when D P >\D Z \. 

Also comparing the experimental specific heat data 
with their theoretical results, Misguich and Sindzingrc 22 
concluded that the measured entropy in the materials 
ZnCu 3 (OH) 6 Cl 2 at T/J = 0.06 is much lower than for 
the pure KLHM, and that the very low temperature spe- 
cific heat may be dominated by impurities. We discuss 
the role of DM interactions in the entropy and specific 
heat of the material. In general, DM interactions should 
lower the entropy as they reduce the manifold of classi- 
cal ground states. We find that the entropy is reduced 
primarily due to D z , whereas D p has a very small effect 
on it. Thus, combining the experimental results on the 
susceptibility and specific heat, we believe that the most 
likely DM parameters for the material are in the range, 



D p / J w 0.2 - 0.3, \D Z \/J « 0.1. However, these results 
may change a little if impurity effects are included in the 
analysis. 

The results of Misguich and Sindzingre^ 2 - further sub- 
stantiate our earlier assertion that the pure KLHM has 
an extended crossover regime^ where the susceptibil- 
ity grows as a power law in inverse temperature and 
the specific heat or entropy is sublinear in temperature. 
From their numerical results, the power law in the sus- 
ceptibility may extend over a full decade in temperature, 
0.1 < T/J < 1. We would like to reiterate that this 
crossover regime is important to address theoretically 
and may well be relevant to the properties of the real 
material. 

This exposition is organized as follows. In Sec. UU we 
discuss the effects of Dzyaloshinsky-Moriya and exchange 
anisotropy on the susceptibility and compare them with 
experiments on ZnCu3(OH)gCl2. In Sec. IIII1 we study 
the consequences of the above mentioned perturbations 
and of quenched dilution on the entropy and specific heat 
and also discuss the constraints on the DM parameters 
for ZnCu3(OH)gCl2 that the specific heat (entropy) mea- 
surements introduce. Finally, the conclusions are pre- 
sented in Sec. HVl 



II. UNIFORM SUSCEPTIBILITY 

In this section, we study the effects of different per- 
turbations on the uniform susceptibility of the spin-i 
Heisenberg model. In a magnetic field h, the Hamilto- 
nian of this model can be written as 



H = J ^2 SjSj + Upert - 9l l B 



(i) 



where J is the exchange coupling, g is the g factor (as- 
sumed to be isotropic), /is is the Bohr magneton, and 
Ti.p er t represents the various perturbations. In what fol- 
lows, we set J = 1 and g/i^ = 1. In the sum, means 
that only nearest neighbor interactions are considered. 
The uniform susceptibility per spin is then given by 



T <9 2 lnZ 



n dhi 



(2) 



h=0 



where Z is the partition function, T the temperature, N 
the number of lattice sites, and a = x,y, z. The molar 
susceptibility Xmoiar, measured experimentally, is related 
to our susceptibility per spin by the relation Xmoiar = 
Cx, where the constant C = N A g 2 ^ 2 B /kJ = 0.3752 g 2 / J 
in cgs units. 

In a kagome lattice (which we assume lies in the x- 
y plane), both out-of-plane (D z ) and in-plane (D p ) DM 
terms are allowed ) 24 ' 25 i 26 ' 27 i 28 



Hdm — 



^(S i xS j ) z + D p .(S i xSj), (3) 
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FIG. 2: (Color online) (a) Orientation of the in-plane and 
out-of-plane DM interactions, [(b)-(d)] Three possible clas- 
sical ground states of the kagome lattice. In each case, the 
system has (b) positive chirality, (c) negative chirality, and 
(d) both positive and negative chiralities, as indicated by the 
signs inside the triangles. 



Since the DM terms break spin rotational symmetry, 
we need to calculate separately the susceptibility with 
field along z (x z ) and in the x-y plane (xp)- The powder 
susceptibility \a is given by Xa = \{^Xp + Xz)- 

With DM interactions, the different S z sectors become 
coupled so that we are able to do NLC calculations only 
up to six triangles^ Unfortunately, the convergence of 
the KLHM with additional DM anisotropy is poor, and 
having only very few terms for the NLC expansion does 
not allow us to perform extrapolations. Hence, we turn 
to ED of clusters with 12 and 15 sites (Fig. [1]), and pe- 
riodic boundary conditions, to study the effects of DM 
interactions. In all our plots, we include both the results 
for 12 and 15 sites. The region where they agree gives 
an idea of the temperature range where finite-size effects 
are small and one can be confident of the ED results. 



A. D z 0, 



D p = 



We first consider the case of a pure out-of-plane DM 
interaction (D z ^ 0, D p = 0). The sign of the D z term 
alternates between the up- and down-pointing triangles 
of the kagome lattice. It can be set by demanding that 
for the up pointing triangle shown in Fig. [2] with corners 
1-2-3, a positive D z multiplies (Si x S2) z , and for the 
down-pointing triangle 3-4-5, a negative D z multiplies 
(S 4 xS 5 ) z . 

The effect of a pure out-of-plane DM interaction is to 
favor the spins to lie in the x-y plane. In that sense, it 
acts like an easy-plane exchange anisotropy. In addition, 
the sign of D z breaks the chiral symmetry of the KLHM. 




FIG. 3: (Color online) Temperature dependence of the (a) in- 
plane and (b) out-of-plane susceptibilities in the presence of a 
pure D z anisotropy. They are compared with \ f° r the pure 
KLHM [computed using the triangle based NLC expansion 
considering up to eight triangles (Ref. |23|) ] and with the ex- 
perimental results (Refs. [7| and |8|) translated to our notation 
(Ref. l29f) . In the inset in (a), we show the powder susceptibil- 
ity and in the inset in (b) the anisotropy. In all plots of DM 
calculations, thick (thin) lines show the ED results of the 15 
(12) site cluster. 



In the classical limit, the ground state of the Heisenberg 
antiferromagnet is highly degenerate. All states satis- 
fying the "120°" rule in each triangle minimize the en- 
ergy. In Figs. [2jb)-[2jd) , we show three possible Neel 
states of the KLHM, where spins lie in the plane and 
which have positive, negative, and mixed chiralities, re- 
spectively. Once the z component of the DM interaction 
is introduced, the degeneracy between the states in Fig. 
[2] is lifted. For D z > 0, the state with negative chiral- 
ity [Fig.[^c)] is favored, while for D z < 0, the state with 
positive chirality [Fig.^b)] is the one that minimizes the 
energy. 

Returning to our spin-half quantum model, in Fig. [3] 
we show the in-plane and out-of-plane susceptibilities as 
a function of temperature for two different strengths of 
D z . (If only D z is present, its sign is irrelevant for 
the thermodynamic quantities we study in this paper.) 
We compare these results with the ones obtained for the 
pure KLHM22 and the ones measured experimentally for 
ZnCu3(OH)eCl2i^ 1 ^ Two features are apparent in these 
plots, which can be related to the increase of in-plane cor- 
relations between spins, as expected from classical con- 
siderations, (i) The D z term suppresses both Xv an d 
Xz with respect to the KLHM result and (ii) Xz becomes 
larger than Xp- Since the experimental result for the pow- 
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der susceptibility is larger than the one for the KLHM, it 
is evident [inset in Fig. 02a)] that a D z term alone can- 
not explain the experiments. The anisotropy produced 
by D z in the susceptibilities is shown in the inset in Fig. 

Hb). 

There is another feature in the plots in Fig.[3]that is no- 
ticeable: finite-size effects set in at higher temperatures 
for Xp than for \z- This confirms that the planar (XY) 
correlations are longer ranged than the ZZ correlations. 



B. 



D p / 0, D z 



We now analyze an in-plane DM interaction D p . ft 
is perpendicular to the bonds and points inward toward 
the center of the triangles. 27 ' 28 This is shown by the ar- 
rows in Fig. EJa). (A different scenario for D p is dis- 
cussed in the Appendix.) D p breaks the rotational sym- 
metry around the c axis, ft also favors classical spin con- 
figurations with a finite z component, producing weak 
ferromagnetismi 26 ' 27 In that sense, D p can act like an 
easy-axis exchange anisotropy. 



X P and Xz with respect to the KLHM result, and that Xz 
becomes larger than Xp- This can be understood consid- 
ering that D p , when trying to produce canting, competes 
with the Heisenberg terms and reduces the in-plane spin- 
spin correlations. 

The enhancement of Xp and Xz seen in Fig. [4] shows 
that a pure D p term can explain the upturn seen exper- 
imentally for the susceptibility. This can be seen in the 
plots of the powder susceptibilities presented in the inset 
in Fig. BJa). Interestingly, the anisotropics produced by 
D z and D p are very similar, as can be concluded by com- 
paring the inset in Fig. BJb) with the one in Fig. OUb). 
They both generate Xz > Xp, and Xz/Xp are °f the same 
order (at least for the intermediate temperatures studied 
here) when D z and D p are of the same order. 

A remarkable difference between the effect of a D z term 
and the effect of a D p term in our finite cluster calcula- 
tion of the susceptibilities is that in the latter, finite-size 
effects are very small as compared to the former one. 
This implies that all correlations are much weaker in the 
presence of a D p term than in the presence of a D z term. 
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FIG. 4: (Color online) Temperature dependence of (a) the in- 
plane and (b) out-of-plane susceptibilities in the presence of a 
pure Dp anisotropy. They are compared with \ f° r the pure 
KLHM (Ref. [23|) and with the experimental results (Refs. Q 
and|f|). In the inset in (a), we show the powder susceptibility, 
and in the inset in (b) the anisotropy. In all plots of DM 
calculations, thick (thin) lines show the ED results of the 15 
(12) site cluster. 

In Fig. 31 we depict the x-y and z susceptibilities as a 
function of the temperature for two different strengths of 
D p £2- As for D z , we compare our results with the ones 
obtained for the pure KLHM 2 ^ and the ones measured 
experimentally. 7,8 We find that a D p term enhances both 



C. Dp 7^ 0, D z 7^ 0, and the experiments 

Once D p 0, in general, thermodynamic quantities 
for D z > start to differ from the ones for D z < 0. 
In Fig. 03 we have plotted the in-plane and out-of-plane 
susceptibilities for a fixed value of D p while increasing 
\D Z \, with D z > [(a) and (b)] and D z < [(c) and (d)]. 
Comparing Fig. 03a) with Fig. 03c) one can see that the 
increase of \D Z \ produces the same effect in Xp no matter 
the sign of D z . For both D z > and D z < 0, the in-plane 
susceptibility is suppressed with respect to its value for 
Dp 0, D z = 0. On the other hand, the effect produced 
by the increase of |-D Z | on Xz is very different depending 
on the sign of D z . Figure G3L) shows that for D z > 0, 
the increase of D z suppresses Xz with respect to its value 
for Dp 0, D z = 0. As seen in Fig. EJd), the opposite 
occurs if D z < 0. The increase of \D Z \ enhances Xz with 
respect to its value for D p ^ 0, D z = 0. 

Given the above results for the in-plane and out-of- 
plane susceptibilities in the presence of D p and D z , one 
can then understand the behavior of the powder averages, 
which is presented in the insets of Figs. 03Ja) and 02b). 
For D p and D z > 0, the increase of D z monotonically 
reduces Xa from its value at D z = 0. On the contrary, 
if Dp 7^ and D z < 0, a small increase of |D Z | enhances 
Xa at lower temperatures, producing a better agreement 
with the experiments than Xa for D p 0, D z = 0. Ulti- 
mately, when \D Z \ ~ D p , the powder average of the sus- 
ceptibility is again suppressed with respect to its value 
for D z = 0. 

Based on the results presented in Fig. 03 we conclude 
that the sharp upturn seen experimentally for the powder 
susceptibilities of the material ZnCu3(OH)gCl2 can be 
understood to be a consequence of DM interactions for 
Dp > \D Z \ and D z < 0. We predict that single crystal 
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measurements should see an upturn in the anisotropy 
when the powder susceptibility departs from the KLHM 
result. Such a behavior has also been seen for spin-| 
kagome system KFe 3 (OH) 6 (S0 4 )2-— 

We should stress that in our theoretical calculations for 
the powder susceptibilities, we have assumed the g fac- 
tor to be isotropic. As depicted in Fig. ODJd) for D z < 0, 
the z susceptibility rises very rapidly, producing a large 
anisotropy Xz/Xp [shown in the inset in Fig. EJd)]. There- 
fore, an expected anisotropic g factor enhanced along z 
will cause an even more rapid rise of Xa than the one 
presented in the inset in Fig.^b), and will lead to agree- 
ment with experiments with a smaller DM anisotropy. 
For example, already for D p — 0.2 J and \D Z \ < D p , the 
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z susceptibility is very similar to the experimental result 
for the powder averages. 

Another factor that could reduce the DM anisotropy 
required to describe the experimental results is the pres- 
ence of a small concentration of impurity spins; 21,22 How- 
ever, as this may vary from sample to sample, we will not 
consider it further in this work. As mentioned in the In- 
troduction, even in the absence of single crystals, recent 
NMR work of Imai et alJ£ may support the relevance of 
DM interaction to ZnCu3(OH)6Cl2. The behavior they 
observe for the main peak and edges of the NMR spec- 
tra is similar to the powder [inset in Fig.[5]Jc)] or Xz [Fig. 
[5jd)] and in-plane [Fig.^c)] susceptibilities, respectively. 
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FIG. 5: (Color online) Temperature dependence of the [(a) and (c)] in-plane and [(b) and (d)] out-of-plane susceptibilities in the 
presence of an in-plane D p = 0.3J anisotropy and different values and signs of the out-of-plane anisotropy. These susceptibilities 
are compared with x f° r the pure KLHM (Ref. l23h and with the experimental results (Refs. and|H). In the inset in (a) and 
(c), we show powder susceptibilities, and in the inset in (b) and (d) the anisotropies. In all plots of DM calculations, thick 
(thin) lines show the ED results of the 15 (12) site cluster. 



D. Exchange anisotropies 

We discuss in what follows the role that exchange 
anisotropies play in the behavior of the magnetic sus- 
ceptibilities. Easy-plane (A > 0) and easy-axis (A < 0) 
exchange anisotropies can be introduced perturbing the 
Heisenberg Hamiltonian with a term 

H EA = Aj2(SfSf + SfSy). (4) 

(i,3) 

Equation (0| breaks the SU(2) symmetry of the 
KLHM, however; it does not couple different S z sec- 
tors. Hence, in the absence of x-y magnetic fields, we 
can perform calculations for out-of-plane susceptibilities, 
entropy, and specific heat using the triangle-based NLC 
expansion summing contributions of up to eight triangles. 



Calculating the in-plane susceptibilities [Eq. @] requires 
introducing x-y magnetic fields, which couples different 
S z sectors. Consequently, NLC calculations for Xv can 
be done only up to six triangles. As we will show later, 
in the presence of exchange anisotropies, NLC results for 
X P are very similar to the ones obtained with ED (15 site 
cluster) down to T ~ 0.35 J. For all other quantities, we 
only present NLC results, which are more accurate. 



1. Easy-plane exchange anisotropy 

Figure[5]shows that an easy-plane exchange anisotropy 
decreases both the in-plane and out-of-plane susceptibil- 
ities with respect to the pure Heisenberg model, and Xz 
becomes larger than Xpi similar to the effect of a pure 
D z term shown in Fig. [3] There are, however, clear dif- 
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FIG. 6: (Color online) Temperature dependence of the (a) 
in-plane and (b) out-of-plane susceptibilities in the presence 
of an easy-plane exchange aniso trop y. Results are compared 
with x for the pure KLHM (Ref. |23T ) and with the experimen- 
tal results (Refs. and|H). In the inset in (a), we show the 
powder susceptibility, and in the inset in (b) the anisotropy. 
In (a) and the insets of (a) and (b), the susceptibilities were 
computed using ED. Thick (thin) lines show the ED results 
of the 15 (12) site cluster. In (b), the z susceptibilities were 
obtained using NLC. Thick (thin) lines are the results of the 
NLC expansion with up to eight (seven) triangles. 



ferences between A > and D z ^ 0. The easy-plane ex- 
change anisotropy produces a large reduction of Xp [Fig. 
EH a)] and Xa [inset in Fig.^a)] at high temperatures, and 
this reduction does not depend strongly on temperature 
down to T ~ 0.3 J. 

The above property of the exchange anisotropy high- 
lights how remarkable DM interactions are. Their effect 
on the susceptibilities is negligible at high temperatures 
(even for large values of D z ) and only onsets as the tem- 
perature is lowered. Another important difference be- 
tween an easy-plane exchange anisotropy and the pure 
out-of-plane DM term is that, as shown in the inset in 
Fig-EDJb), the behavior of the anisotropy Xz/Xp produced 
by the former perturbation is nonmonotonic with tem- 
perature. Hence, susceptibility experiments with single 
crystals should be able to easily distinguish between these 
two types of anisotropics. 



2. Easy-axis exchange anisotropy 

The effect of an easy-axis exchange anisotropy on the 
susceptibility is depicted in Fig. [7] Like D p , A < en- 
hances both Xp an d Xz with respect to x in the pure 
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FIG. 7: (Color online) Temperature dependence of the (a) in- 
plane and (b) out-of-plane susceptibilities in the presence of 
an easy-axis exchange anisotropy. Results are compared with 
X for the pure KLHM (Ref. [231 ) and with the experimental 
results (Refs. and @). In the inset in (a), we show the powder 
susceptibility, and in the inset in (b) the anisotropy. In (a) 
and its inset, we compare results of the NLC expansion with 
up to six triangles (thick lines) with ED of a 15 site cluster 
(thin lines). In (b), the susceptibilities were obtained using 
the NLC expansion with up to eight triangles (thick lines) 
and seven triangles (thin lines). The anisotropies [inset in 
(b)] were obtained using ED of 15 site (thick lines) and 12 
site (thin lines) clusters. 

KLHM case. However, A < has a large effect on Xp 
[Fig. 0a)] and Xa [inset in Fig. Efa)] at high tempera- 
tures. Basically, A < produces a large enhancement of 
the susceptibility that is not strongly dependent on the 
temperature down to T ~ 0.3J. 32 Hence, A < can- 
not provide an explanation for the sharp increase seen 
experimentally in the powder susceptibility. Notice also 
that in presence of an easy-axis exchange anisotropy, the 
in-plane susceptibilities are larger than the out-of-plane 
ones, and their anisotropy [inset in Fig.[7fb)] is nonmono- 
tonic in temperature. 

III. ENTROPY AND SPECIFIC HEAT 

In this section, we to study the entropy (S), 

S=±(lnZ+(H)/n (5) 
and specific heat (C„), 

C v = J= {(H 2 ) (H) 2 ) , (6) 
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of the KLHM in the presence of various perturbations. 
In Eqs. ([5]) and ©, N stands for the number of lattice 
sites and Z for the partition function. 



A. D z / 0, D p = 0, and A > 

We start by considering the effect of an out-of-plane 
DM anisotropy. Results for S and C v and different values 
of D z are shown in Fig. [5] As expected from the quali- 
tative analysis in the previous section, where we argued 
that D z breaks the degeneracy among different possible 
Neel states for KLHM, Fig. [5fa) shows that as the tem- 
perature is reduced, D z reduces the entropy with respect 
to its value for the KLHM. This suppression of the en- 
tropy is accompanied by a large increase in the specific 
heat with respect to the KLHM. While for very small 
values of D z a high-temperature peak can still be seen in 
C v [Fig. Hfb)], by the time D z — 0.3J, any evidence of 
such a peak has disappeared. 
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CO 



0.4 
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of D z [Fig. EJa)], it is only noticeable as the temperature 
is lowered. 

The effect of A > on the specific heat is very dif- 
ferent from that of D z , as seen in Fig. GDJb). Increasing 
A (up to A = 0.3J) only displaces the high-temperature 
peak to higher temperatures almost without changing its 
height. (Eventually, as the exchange anisotropy is further 
increased toward the XY limit, the height of the C v peak 
would decrease.) This means that if one could extract the 
behavior of S and C v from experiments, one could further 
distinguish between DM and exchange anisotropics. 
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FIG. 9: (Color online) (a) Entropy and (b) specific heat as 
a function of the temperature in the presence of an easy- 
plane exchange anisotropy. All results were obtained using 
the NLC triangle-based expansion fRef.l23h considering up to 
eight triangles (thick line) and seven triangles (thin line). 
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FIG. 8: (Color online) (a) Entropy and (b) specific heat as a 
function of the temperature in the presence of an out-of-plane 
DM interaction. For the KLHM, we used the NLC triangle- 
based expansion (Ref. [23h considering up to eight triangles 
(thick line) and seven triangles (thin line). For all plots with 
D z 0, thick (thin) lines show the ED results of the 15 (12) 
site cluster calculation. 

As mentioned in Sec. [Til the easy-plane exchange 
anisotropy [Eq. produces some features that are sim- 
ilar to a pure D z term. In Fig. [§1 we show the effects 
of A > on S and C v . Indeed, like D z , an easy-plane 
anisotropy suppresses the entropy with respect to the 
KLHM. This reduction, however, is apparent at high 
temperatures even if A is small, while in the presence 



B. Dp 0, D z = 0, and A < 

We now turn to the effects of a pure in-plane DM 
anisotropy on the entropy and specific heat. This is de- 
picted in Fig. [TUJ Figures [TOfa) and [TOTb) show that, 
at least down to temperatures ~0.15J, D p has a negligi- 
ble effect on the entropy and specific heat, respectively. 
Considering that the largest D p in Fig. [TO] is 30% of J 
and that such anisotropy produces large changes in the 
uniform susceptibilities (Fig. [4]) , we find this feature re- 
markable. 

We should add that as D p is increased, the small dis- 
placement seen in the high temperature peak [Fig.fTOTb)] 
toward higher temperatures is of the same order as the 
ED finite-size effects. To make that clear, we have also 
plotted Dp = results obtained from the exact diagonal- 
ization of a finite cluster with 15 sites. As better seen in 
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FIG. 10: (Color online) (a) Entropy and (b) specific heat as 
a function of the temperature in the presence of an in-plane 
DM interaction. For the KLHM, we used the NLC triangle- 
based expansion (Ref. [23h considering up to eight triangles 
(thick line) and seven triangles (thin line). For all plots with 
D p 7^ 0, thick (thin) lines show the ED results of the 15 
(12) site cluster calculation. In (b), the extra plot denoted 
by Dp — is the KLHM ED result obtained with a 15 site 
cluster. The inset in (b) magnifies the high-temperature peak 
of C v so that ED finite-size effects become discernible. 



the inset, the ED peak for D p = is slightly displaced to- 
ward higher temperatures than the thermodynamic limit 
result provided by NLC, and that displacement is of the 
same order as the one seen for D p ^ 0. (C v is, in general, 
very sensitive to finite-size effects^) 

The presence of an easy-axis exchange anisotropy has a 
very different effect on S and C v of the KLHM. This can 
be seen by comparing Fig. [TT] with Fig. [TO] The increase 
of |A| increases the entropy, at all temperatures, with 
respect to the KLHM. This is expected since in the Ising 
limit, the system has a finite entropy at zero temperature. 
For C v , what happens is that the high-temperature peak 
moves to lower temperatures. The height of the peak 
almost does not change up to A = — 0.3J, although it 
eventually decreases as A approaches -122- 



C. Dp^O, D z ^0 

If both DM terms are present in the system, we find 
that the deviations of S and C v from the KLHM re- 
sult (in the range of temperatures discussed in this 
work) are mainly determined by the value of D z , al- 
most independent of the value (and sign) of D p (up to 
Dp ~ 0.3J). Hence, at least at intermediate and high 
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FIG. 11: (Color online) (a) Entropy and (b) specific heat as 
a function of the temperature in the presence of an easy-axis 
exchange anisotropy. All results were obtained using the NLC 
triangle-based expansion f Ref. |23|) considering up to eight tri- 
angles (thick line) and seven triangles (thin line). 



temperatures, S and C v are quite insensitive to the ex- 
istence of an in-plane DM interaction. This is further 
discussed in a later section in comparison to experiments 
on ZnCu 3 (OH) 6 Cl 2 . 



D. Quenched dilution 

As discussed in Ref. [ill , another important perturba- 
tion of the KLHM is the presence of quenched dilution. 
Quenched dilution could be generated in the materials 
ZnCu3(OH) 6 Cl2 due to the substitution of Cu sites in the 
kagome planes by Zn. The missing spins on the lattice 
could create local moments in the singlet background and 
cause a Curie-like susceptibility to arise as the temper- 
ature is lowered. However, we have shown that at least 
down to temperatures T ~ 0.3 J, the only effect that such 
a dilution has on the susceptibility is to suppress it with 
respect to the KLHM result^ In this sense, quenched 
dilution has the opposite effect of impurity spins as the 
latter enhance the susceptibility. 

Here, we present studies of the effects of quenched dilu- 
tion on entropy and specific heat of KLHM. Our calcula- 
tions are performed using the triangle-based NLC expan- 
sion considering up to eight triangles. If c is the dilution, 
we assume that at each site we have a hole with proba- 
bility c and a spin with probability 1 — c. The holes are 
fixed in their position and extensive quantities are aver- 
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aged over all possible configurations C using the relation 

(o) = J2p(c)o(C), (7) 

c 

where 

P(C) = cF h (l-c) N ' (8) 

is the probability of the configuration C with Nh holes 
and iV s spins. 
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FIG. 12: (Color online) (a) Entropy and (b) specific heat as 
a function of the temperature in the presence of quenched 
dilution. All results were obtained using the NLC triangle- 
based expansion (Ref. |23| ) considering up to eight triangles 
(thick line) and seven triangles (thin line). 

The entropy and specific heat for several values of hole 
concentration are shown in Fig. [T2J We note that at 
these intermediate temperatures [Fig. [T27 a)], holes sim- 
ply lower the entropy at all temperatures. In the case of 
the specific heat, they just displace the high temperature 
peak to lower temperatures, almost without changing its 
height. 



E. Entropy difference from kagome lattice 
Heisenberg model and implication for 
ZnCu 3 (OH) 6 Cl 2 

To conclude the section on specific heat and entropy, 
we study the entropy difference between the KLHM and 
the model with different DM parameters. Shown in Fig. 
fl~S1 is AS*, given by 

AS = S{D P = 0,D Z = 0)- S(D P , D z ). 



It is the reduction in entropy due to DM interactions. In 
Fig. 021 one notices that the entropy reduction is deter- 
mined primarily by D z and D p plays a small role. 

The diamond in Fig. [13] represents the minimum dis- 
crepancy between the pure KLHM and ZnCu3(OH)6Cl2 
at T I J — 0.06 as determined by Misguich and Sindzingre 
(MS) in Ref. The experimental data at higher tem- 
peratures are likely dominated by phonons. From Fig. 
[TBI we conclude that \D Z \/J is likely to be about 0.1 in 
ZnCu3(OH)6Cl2. Hence, while the sharp increase in the 
susceptibility discussed in Sec. [TT] allows us to make an 
estimate of the possible values of D p (D p / J w 0.2 — 0.3), 
the entropy reduction allows us to get an estimate of \D Z | . 




FIG. 13: (Color online) Difference between the entropy of the 
pure KLHM and the entropy of the KLHM in the presence 
of DM anisotropies. Thick dashed lines show the ED results 
of the 15 site cluster and the thin black lines result of 12 site 
cluster. The diamond at T — 0.06J depicts the minimum dif- 
ference between the pure KLHM and the experimental result 
as obtained by MS in Ref. [H. 



IV. CONCLUSIONS 

In this work, we have studied in detail the effect of 
perturbations such as in-plane and out-of-plane DM in- 
teractions, exchange anisotropies, and quenched dilution 
on the KLHM. We have focused here on the effect such 
perturbations have on magnetic susceptibilities, entropy, 
and specific heat. We first summarize our theoretical 
findings: 

(i) In the presence of a pure out-of-plane DM term 
(D z ^ 0, D p = 0) both in-plane (x P ) and out-of-plane 
(Xz) susceptibilities are suppressed with respect to the 
KLHM result, and \z becomes larger than \p- However, 
the susceptibility anisotropy Xz/Xp is n °t large (< 1.1) 
down to T ~ 0.5J for D z < 0.3J. On the other hand, 
a D z term suppresses the entropy as the temperature is 
lowered and produces an increase of the specific heat at 
intermediate temperatures. 

(ii) In the presence of a pure in-plane DM term (D z = 
0, D p ^ 0), both in-plane (x P ) and out-of-plane (xz) 
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susceptibilities are enhanced with respect to the KLHM 
result and, as for pure D z , \z becomes larger than Xp- 
The susceptibility anisotropy is not large (< 1.1) down 
to T ~ 0.25J for D p < 0.3 J. A D p term (< 0.3 J) has a 
negligible influence on entropy and specific heat down to 
T ~ 0.15J. 

(iii) When both D p and D z are present, the suscep- 
tibility becomes sensitive to the sign of D z . (A) For a 
constant value of D p , the increase of D z for D z > sup- 
presses both the in-plane and out-of-plane susceptibili- 
ties with respect to the D z = value. The susceptibility 
anisotropy is of the same order as when only D p or D z 
is present. (B) Also keeping D p constant, the increase 
of \D Z \ for D z < suppresses the in-plane susceptibility 
but enhances the out-of-plane susceptibility. 

(iv) In the presence of both D p and D z , the entropy 
and specific heat are mainly determined by the modulus 
of D z . D p and the sign of D z play a relatively small role. 

(v) The presence of an easy-plane exchange anisotropy 
suppresses both Xp an d Xz from their KLHM values, and 
Xz becomes larger than \v Such anisotropy has a large 
effect even at high temperatures and does not produce 
deviations from the KLHM that are strongly tempera- 
ture dependent. In addition, it generates susceptibility 
anisotropies that are nonmonotonic and weakly depen- 
dent on temperature down to T ~ 0.25J. A > also 
suppresses the entropy with respect to the KLHM. The 
high-temperature peak of the specific heat is slightly dis- 
placed toward higher temperatures almost without mod- 
ifying its height (for A < 0.3 J). 

(vi) The presence of an easy-axis exchange anisotropy 
enhances both Xp an d Xz from their KLHM values, and 
Xp becomes larger than Xz- A < has a large effect 
on the high-temperature values of Xp an( i a ^ so does not 
produce deviations from the KLHM that are strongly 
temperature dependent. It also leads to susceptibility 
anisotropies that are nonmonotonic and weakly depen- 
dent on the temperature down to T ~ 0.25J. A < 
also enhances the entropy with respect to the KLHM 
result. The high-temperature peak of the specific heat 
is slightly displaced toward lower temperatures without 
much change in height (for |A| < 0.3 J). 

(vii) At intermediate and high temperatures (T > 
0.3J), quenched dilution has been shown to suppress the 
uniform susceptibility with respect to the KLHMf^ We 
have discussed here that it also reduces the entropy for 
all temperatures T > 0.3J. In the case of the specific 
heat, the effect of dilution is to displace the high temper- 
ature peak toward lower temperatures without affecting 
its height (at least when c < 0.2). 

We now discuss our conclusions with regard to the ob- 
served properties of the material ZnCu3(OH) 6 Cl2: 

(i) The observed susceptibility shows large enhance- 
ment with respect to the KLHM, which has a sudden 
onset below T = J/2. This kind of behavior is only 
compatible, within the models studied, with a large 
D p w 0.2 - 0.3 and &D Z <0, with \D Z \ < D p . 

(ii) Misguich and Sindzingre^ have previously con- 



cluded that the experiments on ZnCu3(OH)gCl2 show 
that there is a large reduction in entropy with respect 
to KLHM at T/J = 0.06. This reduction is at least 
0.05 and may be larger when impurities and phonons are 
taken into account. The comparison with DM anisotropy 
calculations shows that this implies \D Z \/J w 0.1. 

Based on these results, our overall conclusion for the 
parameters of the material ZnCu3(OH)gCl2 is D p /J in 
the range 0.20 - 0.30, \D Z \/J w 0.1, and J « 170K. 
All these numbers could change if there are substantial 
impurity contributions present. Future experiments on 
the effects of anisotropy can resolve these issues. 

On the theoretical side, an important question that re- 
mains open is the nature of the low-temperature phase(s) 
of the KLHM in the presence of DM terms. In clas- 
sical systems, it has been shown that ground state is 
ordered ) 26 ' 27 and there is a finite temperature phase tran- 
sition in which the critical temperature depends on the 
values of D z and D v *& However, one should keep in 
mind that for classical systems, even in the absence of 
DM anisotropy, the system orders as T goes to zeroi 34 ' 35 
Quantum effects for the spin-j case, and their relation to 
the experimental absence of any order down to 50 mK, 
still need to be elucidated. 
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Note added on proof. -It has been argued in recent 
experiment s 36 ' 37 ' 38 that contrary to the original expec- 
tation a rather large concentration of intersite mixing 
(Cu/Zn) impurities [c s ~ 6 — 10%) may be present in 
the Hcrbertsmithite material. One natural question that 
has not been addressed in these works is whether such 
a large impurity concentration can be accommodated to 
reproduce the intermediate and high temperature mag- 
netic susceptibility of the KLHM. In Ref. [ll| we have 
shown that subtracting the contribution of c s = 4.5% 
free impurity spins (J — 200 K) to the experimental re- 
sults of Refs. and [H one can reproduce the magnetic 
susceptibility of the KLHM down to T = 0.3J (The 
free impurity spin concentration reported in Ref. [ll| was 
incorrect by a factor 3/2, which means that in Fig. 2 of 
that reference one should read c = 0.045 and c = 0.09 in- 
stead of c = 0.03 and c = 0.06, respectively.) However, a 
sharp rise in the susceptibility remained below T = 0.3J 
that was not expected in the KLHM. It was shown later 
by Misguich and Sindzingre 2 ^ that adding a small fer- 
romagnetic coupling between impurities (c s = 3.7% and 
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J = 190 K) allows one to reproduce the magnetic suscep- 
tibility of the KLHM down to T = 0.2 J, with \ starting 
to decrease when T ~ 0.1J. 

From the theoretical results in Refs. [HI and one 
could conclude that after subtracting the large contribu- 
tion of c s ~ 6 — 10% impurity spins (Cu) to the experi- 
mental results one would obtain a magnetic susceptibility 
that is incompatible with the KLHM at intermediate and 
high temperatures. However, such a large concentration 
of impurity spins can be accommodated if one consid- 
ers the effect of the non-magnetic impurity counterpart 
(Zn) that is present in the kagome planes. As shown 
in Ref. [ll] such impurities reduce the magnetic suscepti- 
bility with respect to the KLHM. Figure [14] shows that 
after subtracting a c s ~ 6% free impurity spin contribu- 
tion from the experimental results one can reproduce the 
susceptibility of the KLHM with a c = 6% of quenched 
nonmagnetic impurities. 
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FIG. 14: (Color online) NLC results for the KLHM with 
(c = 6%) and without (c = 0) nonmagnetic impurities are 
compared with the experimental results after a c 3 — 5.5% 
and c s = contribution of free impurity spins is subtracted, 
i.e., in the latter case we have plotted Xmoiar/C — c s /(4T) 
(where C was defined in Sec. [TTJ . In the presence of impuri- 
ties we obtain J — 210 K, g = 2.37 as opposed to J = 170 K, 
g = 2.33 in their absence. NLC results for the eight(seven) 
triangle based expansion are plotted as thick(thin) lines. 



The results presented in Fig. [LA show that a large con- 
centration of intersite mixing (Cu/Zn) impurities is com- 
patible with the KLHM at intermediate and high tem- 
peratures. However, whether impurities are the main 
contribution to the magnetic susceptibility at interme- 
diate and low temperatures, as opposed to DM interac- 
tions, still needs to be clarified measuring the anisotropy 
in the in-plane and out-of-plane susceptibilities for tem- 
peratures where x clearly departs from the pure KLHM 
result. 



0.15 
0.13 



*f 0.11 



(a) 



KLHM 

D =0 




FIG. 15: (Color online) Temperature dependence of the (a) 
in-plane and (b) out-of-plane susceptibilities in the presence 
of a D p = 0.3 J anisotropy (with D p alternating from triangle 
to triangle as explained in the text) and different values of D z . 
The results are compared with x f° r the pure KLHM (Ref. 
I23T ) In the inset in (b), we show the anisotropy produced by 
these DM terms. In all plots of DM calculations, thick (thin) 
lines show the ED results of the 15 (12) site cluster. 



APPENDIX A: DZYALOSHINSKY-MORIYA 
INTERACTIONS UNDER A DIFFERENT 
LATTICE SYMMETRY 

In Sec. UH we discuss the effects of the DM anisotropy 
allowed by the symmetry of ZnCu 3 (OH) 6 Cl2 on the sus- 
ceptibility of the KLHM. In ZnCu 3 (OH) 6 Cl 2 , the bonds 
between Cu 2+ ions (which have an oxygen atom in the 
middle) are distorted away from the kagome planes, and 
the direction of this distortion alternates from triangle 
to triangle. This symmetry of ZnCu3(OH)eCl2 sets the 
direction of D p , which within our notation is to be al- 
ways pointing toward the center of the triangles i 26 i 27 [In 
Fig. O the Dp that multiplies (Si x 82)^ is (0,-D p ,0), and 
the one that multiplies (S 4 x S 5 ) y is (0,— D p ,0).] If the 
bonds between the Cu sites would not be distorted at all, 
a perfect kagome lattice would be embedded in three di- 
mensions, because of the symmetry of the lattice D p = 0, 
and only D z could be different from zero i 24 i 25 

For completeness, we briefly discuss here how the sus- 
ceptibility of the KLHM would behave in a material 
where the bonds between magnetic ions are all distorted 
in the same direction away from the kagome planes. In 
this case, the D z terms in Eq. ([3]) are identical to the 
ones we have considered for ZnCu 3 (OH) 6 Cl2, but the D p 
vectors will, within our notation, alternate pointing in- 
ward and outward of the up-pointing and down-pointing 
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triangles, respectively. That scenario seems to be very 
interesting within the discussion in Ref. 

In Fig. HSl we depict the x-y and z susceptibilities as 
a function of the temperature for D p = 0.3 J and three 
different values of D z . One can see in Fig. [15] that, for 
D z = 0, the kind of D p anisotropy discussed in this ap- 
pendix has almost no effect on the susceptibility of the 
KLHM, at least for the temperatures considered in this 
work. The inset in (b) shows that it also does not gen- 
erate any asymmetry between \p an d \z- Introducing 
a finite D z , either positive or negative, only suppresses 



(in a very similar way independent of the sign) both in- 
plane and out-of-plane susceptibilities with respect to the 
KLHM. (Similar to the discussion in Sec. IIIIl the D p 
anisotropy discussed in this appendix is almost irrelevant 
to the entropy and specific heat.) 

Hence, in kagome lattice materials with a crystal 
symmetry different from that of ZnCu3(OH)gCl2, D p 
anisotropics will produce a very different behavior of the 
susceptibility, which will not be enhanced with respect 
to the one of the KLHM. 
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